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ABSTRACT 
This  theoretical  study  of  the  trajectories  of  charged  particles 
supplements  an  experimental  project  of  molecular-ionic  rearrangements 
using  the  magnetic  focusing  properties  of  a  beta-ray  spectrometer.    The 
experimentally  measured  magnetic  field  was  analytically  represented 
by  a  twelfth-order  polynomial.    This  field  is  axially  symmetric ,  but 
non-homogeneous  otherwise.    The  particular  particle  of  concern  was 
the  hydrogen  ion,  H2+.    The  trajectory  of  this  particle  was  computed 
from  a  second-order  differential  equation,  assuming  values  for  the 
kinetic  energy  and  initial  angle  of  the  particle,  and  the  magnetizing 
current.    The  solution  was  obtained  by  numerical  integration  using  a 
CDC  1604  digital  computer.    The  distinctive  feature  of  these  calcu- 
lations, in  contrast  to  those  used  normally  for  a  beta-ray  spectrometer, 
is  the  large  scattering  (initial)  angle,  about  45°. 
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TABLE  OF  SYMBOLS 

A  The  vector  potential 

B  The  magnetic  flux  density  of  the  field 

B(i)  The  coefficients  computed  for  the  polynomial  which  describes 
the  magnetic  field  on-axis 

E  The  kinetic  energy  of  the  rebound  particle,  a  function  of  the 
scattering  angle 

E0  The  energy  of  the  incident  particle 

Fi,  Gi  Functions  of    z    which  are  described  later  in  detail 

H0(z)  The  axial  component  of  the  field  on  axis 

H1  The  i-th  derivative  of  the  function  H0(z) 

IB  The  magnetization  current 

k  =  p/q 

L  The  intercept  distance  measured  on  the  z-axis 

0  The  azimuthal  angle  in  cylindrical  coordinates 

p  The  momentum  of  the  rebound  particle  (H2+) 

q  The  magnitude  of  the  electronic  charge 

r  The  radial  distance  in  cylindrical  coordinates 

6  The  scattering  angle  measured  in  the  laboratory  frame 

z  The  axial  distance  in  cylindrical  coordinates 


1  .      Introduction 

This  presentation  is  part  of  an  experimental  project  presently  being 
conducted  at  the  U.S.  Naval  Postgraduate  School  to  verify  a  classical 
theory  developed  by  Bates,  Cook,  and  Smith  [l  ]  of  the  capture  of  a 
light  ion  or  atom  from  a  target  system  by  a  fast  projectile. 

Specifically,  the  theory  involves  a  reaction  of  the  type: 

H+  +  CH4       -      H+H  +  CH3 

and  is  essentially  a  sequence  of  binary  collisions.    The  process 
ensues  with  the  bombardment  of  methane  gas  by  low  energy  protons, 
wheeorby  the  resulting  rebound  proton  and  escaping  hydrogen  ion  join  to 
form  the  product  H2+  .     The  angle  of  departure  (or  scatter)  is  shown  to 
be  in  the  vicinity  of  45  degrees,  and  is  labeled  Q. 

Since  the  expected  cross  section  will  be  quite  low  (^10-18  cm2), 
a  detector  could  not  be  utilized,  as  is  standard  procedure,  by  placement 
at  the  scattering  angle.     [5  ] 

Therefore,  the  beta-ray  spectrometer  was  chosen  due  to  its  prop- 
erty which  will  cause  all  H2+  ions  of  the  same  departure  angle  and 
energy  to  converge.    This  study  determines  that  point  of  convergence, 
or  focus . 
2.      Discussion 

It  is  known  from  elementary  electromegnetic  theory  that  a  magnetic 
force,  q{v  x  B},  is  produced  on  a  charged  particle  moving  at  an  angle 
to  the  direction  of  the  field.    This  force  causes  the  particle  to  sustain 
a  change  in  velocity  (direction  only).    The  supposition  is,  for  this 
field,  that  a  particle  which  begins  its  motion  on  the  axis  of  the  field 
and  at  an  angle  with  that  field  axis ,  will  travel  a  curvilinear  path  and 
ultimately  intersect  the  axis  of  the  field  for  appropriate  values  of  the 
field  strength,  energy  and  initial  angle  of  the  particle.    This  can  be 
readily  shown  for  a  homogeneous  field. 

The  point  at  which  the  particle  returns  to  the  axis  is  called  the 
intercept  point,  and  the  distance  from  the  origin  of  the  trajectory  to 
the  intercept  point  is  the  intercept  distance  (L) .    This  is  where  the 


detector  can  be  placed  to  analyze  the  particulate  beam .    The  experimental 
usefulness  is  that  if  the  intercept  distance  is  known  along  with  the 
field  strength,  and  the  energy  of  beam  E0/  then  the  initial  angle  may  be 
evaluated. 

The  theoretical  problem,  which  is  herein  shown,  is  the  evaluation 
of  the  intercept  distance  assuming  values  for  the  scattering  angle. 

Thus,  the  following  intermediate  problems  arose: 

A.  A  means  to  express  the  field,  or  more  precisely,  a  continuous 
function  which  would  represent  the  field  for  all  points  in  the  area  of 
interest. 

B.  A  set  of  equations  to  specify  the  motion  of  the  particle  to 
incorporate  the  energy  and  charge  of  the  particle,  as  well  as  the  inter- 
action of  the  magnetic  field. 

C.  A  numerical  method  to  solve  a  differential  equation  which 
would  result  from  the  equations  in  (B)  above. 

D.  Interpretation  and  extrapolation  from  sets  of  trajectories  to 
obtain  these  derivatives: 

(-^-)  —  The  change  in  intercept  distance  with 

\  dE/0,  IB  respect  to  changes  in  energy 

/aiA 

\b6JlQ  ,  E0     —  The  change  in  intercept  distance  with 
respect  to  changes  in  angle 


\MB/  E0. 


The  change  in  intercept  distance  with 
respect  to  changes  in  magnetization  current 


3.      The  Magnet  and  its  Field 

The  magnet  which  will  be  used  for  this  project  was  previously  in 
the  USNPGS  beta-ray  spectrometer.    It  is  essentially  composed  of  a 
normal  section  of  a    right  circular  cylinder,  and  wound  with  cylindrically 
concentric  layers  of  copper  wire.    The  outer  disks  and  cylindrical 
center  section  which  form  the  container  for  the  layers  of  wire  are 
copper  as  are  the  cooling  water  tubes.    The  figure  (Fig.  1)  on  the 
next  page  shows  the  magnet  with  its  dimensions . 
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Fig  .  1  .     The  Magnet 
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The  experimental  field  strength  of  this  magnet  on-axis,  H0,  appears 
as  a  Gaussian  curve  as  shown  in  Fig.  2.    The  quantities  in  Appendix  1 
for  the  values  of  the  field  are  shown  in  the  positive    z    direction.    The 
maximum  intensity  of  the  field  on-axis  is  at  the  geometrical  center  of 
the  magnet,  at  the  point  formed  by  the  intersection  of  lines  ii  and  jj . 
B  has  cylindrical  symmetry  about  the  line    jj ,  and  is  symmetrical  about 
the  center  plane. 

In  order  to  more  easily  discuss  the  field,  the  system  of  cylindrical 
coordinates  should  be  established  at  this  point.    The  origin  is  defined 
to  be  at  the  intersection  of  the  lines    ii    and    jj .    The  line    jj    will  be 
the    z    axis,  and  the  radial  distance,  r,  will  be  measured  perpendicularly 
from  it.    The  angular  component,  0,  is  as  in  (r,  0,  z) ,  a  right  hand 
coordinate  system. 

Measurements  were  made  on  this  field  for  sarious  values  of  field 
current  by  P.  J.  Kelly  in  May  1965  with  an  axial  Hall  Probe,  and  this 
data  has  been  compiled  in  Appendix  1 .    The  axial  component  of  the 
field  intensity  was  measured  initially  on  the  primary  axis  of  the  magnet 
(z  ax£s)  for  a  distance  of  66  cm.  on  either  side  of  the  origin  at  4  cm. 
intervals.    This  was  done  for  various  current  settings  between  2  and 
14  amperes . 

Also,  the  axial  component  of  the  field  was  measured  off  axis  at 
four  values  of    z,  and  for  8  values  of    0,  for  each  value  of    z.    The 
radius    r   was  varied  from  zero  to  seven  cm. 

The  conclusions  concerning  the  field  were  that: 

(1)  The  field  strength  is  linear  with  current  over  the  entire 
range  of  current  settings  within  an  accuracy  of  ±0.1%. 

(2)  The  field  is  symmetrical  on  either  side  of  the    z  =  0 
plane  to  ±0.01%. 

(3)  The  field  is  azimuthally  symmetric  to  ±0.01%. 

To  find  the  analytic  expression  for  the  field  on  axis,  H0  (z) ,  where 
z    assumes  positive  values  only,  a  polynomial  was  fit  to  the  data  by 
the  method  of  least-squares.    This  data  is  contained  in  Appendix  1  .    A 
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Fig.  2.    The  Magnetic  Field 
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twelfth  order  polynomial  was  chosen,  and  provided  sufficient  precision 
in  the  form: 

H0  (z)  =  Bj_  +  B2  z  +  B3  z2  +  B4  z3  +  B5  z4  +  .    .    .  +  Bl3  z12     .  (3- 1 ) 

The  thirteen  coefficients  were  determined  by  the  computer  and  a  USNPGS 
subroutine  LSQPOL,  which  is  in  PROGRAM  CURVE  included  in  Appendix 
6. 

At  this  point ,  a  comparison  was  made  of  the  experimental  field 
values  and  those  values  which  would  result  from  the  polynomial  field 
expression. 

To  make  this  eheck  of  the  values  of  the  field  off-axis,  the  following 
steps  were  taken: 

(1)  The  vector  potential  was  computed  for  all  values  of    r    and    z 
using  the  derived  polynomial. 

(2)  The  axial  component  bf  the  field  was  derived,  using  the 
vector  potential,  and  evaluated  for  various  values  of    r    and    z. 

This  computation  is  shown  in  detail  in  Appendix  1 .    The  agreement  be- 
tween experimental  and  theoretical  values  was  again  quite  good. 
4.      Theory  of  the  Trajectory 

The  H3+  ion  has  a  positive  charge  eq<yal  in  magnitude  to  that  of  the 
electron,  and  mass  essentially  equal  to  twice  the  mass  of  a  single 
proton.    The  mass  of  the  single  electron  is  ignored. 

The  path  of  this  particle  can  be  represented  by  this  equation: 


hfHtk,-A,]-r,w[t]  +A[lr 


(4-1) 


where  A    is  the  vector  potential,  and    k    is  defined  as  (mv/q) .    The 
prime  above  indicate  differentiation  with  respect  to    z .    A  discussion 
of   A,  its  derivatives,  and    k    will  follow. 

This  equation  for  any  charged  particle  can  be  derived  from  the 
equations  of  motion  for  that  particle  which  begins  its  path  on  the 
primary  axis  of  the  magnetic  field.    This  equation  is  derived  in  Appendix 
4. 
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The  energy  of  the  H2+  particle  after  scattering  is    E,  and  can  be 
shown  to  be: 

E  =  2  E0    cos2    9  (4-2) 

where    E0    is  the  kinetic  energy  of  the  particle  which  is  incident  on  the 
target  methane.    The  angle    6  is  again  the  angle  of  scatter. 

A  description  of  the  collision  process  is  as  follows: 
The  incoming  proton  strikes  one  of  the  hydrogens  in  the  methane  molecule, 
This  is  the  first  of  the  two  binary  collisions,  and  the  only  one  which  is 
assumed  to  be  totally  elastic.    The  initial  particle  moves  off  at  an 
angle,  called  the  departure  angle;  and  the  second  particle  (essentially 
a  proton  and  an  electron)  moves  towards  the  carbon  atom  at  an  angle 
equal  to  90°  -  6  .     Since  the  carbon  atom  is  assumed  to  be  infinitely 
massive,  compared  to  the  hydrogen  atom  approaching  it,  this  second 
particle  is  redirected  in  the  second  binary  collision  in  the  direction  of 
the  recoiling  proton . 

A  schematic  of  the  collision  process  is  as  follows: 


@ 


In  the  laboratory  frame,  the  initial  energy  of  the  system  is: 

E0  =  ^m1  (vj2  (4-3) 

and  the  final  energy  is: 

E  =  £(2m)(v')2  (4-4) 

Conservation  of  momentum  yields,  for  the  first  collision: 

Px      =  mv    =    mv1  cos  0  +    mv"  sin  8  (4-5) 
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and, 

p  y       =  mv'  s  in  6  -  mv"  cos  6=0  (4-6) 

The  last  equation  results  in 

v"  =  v'  tan  6  (4-7) 

and  substituting  into  the  equation  (4^5) 

"v1  =  v"cos  e  .  (4-8) 

With  this  relationship  between  the  velocities,  the  final  energy  after 
the  two  collisions  becomes 

E  =  i(2m)  (v1)2  (4-4) 

=  |  (2m)  v2  cos2  6 

=  2E0  cos2  6 
E  =  2E0  cos2  6  (4-2) 

If  both  collisions  are  elastic  in  nature,  then  it  is  quite  easy  to  show 
that    0  must  be  45°  .    The  y-component  of  momentum  produces  the 
equation 

Py       =  mv"'  sin  0-  mv1  cos  0=0  (4-9) 

at  the  first  collision.    And  it  can  only  be  satisfied  when 

sin  6=  cos  0 
or,  when  0=  45°. 
5.      Expression  for  the  Vector  Potential 

The  fact  that  the  field  is  cylindrically  symmetric  carries  with  it 
the  implication  that  there  is  no  azimuthal  component,  since  the  curl 
B  =  zero.  This  field  can  be  described  by  a  vector  potential  which  is 
always  in  the  azimuthal  direction. 

Using  the  calculated  expression  for  magnetic  field  on-axis,  H0(z), 

the  vector  potential  expansion 

A(r,z)  =-L_  H0(z)  --^(-^)2  H0(z)  +-jL   (  -M*  H0(z) 

2  16      Sz  384       dz  (5-1) 

2n 


n>n    f^\2n+y_A_yn  how 

n!    (n+1)!       V  2/         \bz  J 
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is  shown  to  be  for  any  point  (r,z): 

A=  i_[G1]-^l[G3]+-ll[G3]--^[G4]+    £ [G5] 

2  16  384  18432  1474560 

-       ^ [G6  ]  +     £1 [G7  ] 

176947200  29727129600 

where    I        ]        H0(z)    is  the  nth  partial  derivative  of  the  function 
v  bz) 

H0(z)  with  respect  to  z,  and  where  Gx  =  H0(z),  and  G2  =[-^-|     H0(z),  etc. 

\bzj 

6.    The  Computer  Solution 

It  remains,  at  this  point,  to  find  the  solution  of  the  differential 
equation.    The  method  chosen  for  this  task  was  a  process  of  numerical 
integration  employing  a  digital  computer.    The  subroutine  INTEG-1  of  the 
USNPGS  Computer  Library  was  selected  because  of  its  adaptability  to 
this  particular  problem.    The  program  in  which  INTEG-1  is  contained 
is  called  PROGRAM  ORBITS,  and  is  listed  as  Program  C  in  Appendix  6. 

The  first  four  statements  of  this  program  (ORBITS)  are  standard 
according  to  the  instructions  for  use  of  INTEG-1  .    The  three  comment 
cards  that  follow  describe  ENER,  OKAY,  and  P,  which  are  names  given 
various  quantities  used  in  the  solution. 

B(l)  through  B(13)  are  the  coefficients  used  in  the  polynomial 
as  obtained  by  PROGRAM  CURVE.  The  functions  F(l)  through  F(6), 
and  G(l)  through  G(7)  are  explained  in  Appendix  5. 

VECDZ  is  the  symbol  for  the  partial  derivative  of  the  vector 
potential  with  respect  to    z.    VECDR  represents  the  partial  derivative 
of  the  vector  potential  with  respect  to    r.    VECPT  is  the  vector 
potential. 

A  second  order  differential  equation  can  be  reduced  to  a  system 
of  two  first  differential  equations.    Applying  this  technique  to 
equation  (4-1)  yields 

[r,  j  .  ,   [»_^l  +  (r,,3  [*  jj-J.  (r,^|f]  .[*  -jf  3 

[k2  -  A2]       [k2  -  A2]      [k2  -  A2]    [k2  -  A2] 

(6-1) 
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[r]"=r'  (6-2) 

Letting 

x(l)  ■  r  and    x(2)  =r'   , 

the  equations  (6-1)  and  (6-2)  become 

[x(2)]'  =  [x(2)]FUNC4  +  [x(2)]3FUNC4  -  [x(2)]2FUNC5  -  FUNC5 

Cx(l)]'  =x(2) 
where 

FUNC4 


LA 

BA 
Bz 

] 

[k2 

-A3 

] 

[A    - 

BA 

Br 

3. 

FUNC5  = 

[k2  -  A2  ] 

Further,  FUNC4  and  FUNC5  can  be  written  in  terms  of  FUNC1 , 
FUNC2  and  FUNC3  as  follows: 

FUNC1  =  k8  -  A2 

FUNC2  =    A  liL 

Bz 

FUNC3  =  A  -M- 
Br 

and 

FUNC4=    FUNC2 
FUNC1 

FUNC5  =    FUNC3 
FUNC1 

OKAY  is  obtained  in  the  following  way.    This  is  the  machine 
language  expression  for    k2    where  the  energy  can  be  entered  in  eV.    In 
the  MKS  system,  k  must  have  the  unit  of  webers/meter.    This  is  the 
same  as  A,  since  B  has  the  units  of  webers/meter2,  and  B  =  7x  A.- 

k  =    m  v      in  units  of    (kilogram) (meter/sec)     or    webers/meter. 
q  (coulombs) 


It  would  be  desirable  to  have    k    expressed  in  terms  of  the  energy 
of  the  particle,  E,  and  in  the  units  of  eV,  as  a  great  savings  in  con- 
version would  be  realized. 

Using  the  non-relativistic  expression  for  the  energy  E,  the 
following  is  obtained: 

E  =_e!_  ;.  p  =  (2mE)£ 

2m 

But,  since  k2  is  needed 

k2  =     P2       =  JmE_    =    2m£'  if  E'  is  expressed  in  eV. 

q2  Q2  q 

And,  k2  becomes 

['  2 (mass  of  Hs+ )    1 
charge  of  electron! 

Inserting  the  values  of  the  constants,  k2  becomes  the  quantity  OKAY: 


I      1.6021  x  1CT19  J 
=    {4.17576  x  10"8  }  E' 


or 

OKAY  =  (4.17576E  -  08)*(ENER),  as  expressed  in  the  program. 

The  method  of  solution  employed  by  INTEG-1  is  a  fourth  order 
Runge-Kutta  numerical  integration  approximation  of  the  ordinary  differ- 
ential equations . 

The  computation  is  performed  in  the  MKS  system  of  units. 
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TO  USE  PROGRAM  ORBITS  (the  main  trajectory  program) 

This  program  is  written  in  the  language  of  FORTRAN- 60.    Its 
successful  employment  depends  on  the  correct  use  of  the  data  cards, 
which  are  the  last  eight  cards  of  the  program  itself.    These  eight  cards 
contain  the  input  information,  and  the  quantities  desired  as  output. 
An  example  of  the  data  cards  may  be  seen  on  the  last  page  of  Program 
C  in  Appendix  6. 

Data  card  #1 :    This  card  is  used  for  identification,  and  columns  2-32 
are  available  for  this  use.    This  title  appears  on  all  data  output,  as 
well  as  graphs . 

Data  card  #2:    This  card  contains  the  number  of  "runs",  or  trajectories, 
to  be  processed  with  a  maximum  of  9.    This  number  is  entered  in 
column  1  of  the  card.    Note,  however,  that  only  one  graph  is  possible 
for  the  entire  program  . 

Data  card  #3:  The  number  of  coupled  first-order  differential  equations 
is  specified  by  this  card  (maximum  of  30).  In  this  case,  there  are  two 
and  this  number  is  posted,  right  justified,  in  columns  1  and  2,  as  02. 
Data  card  #4:  The  initial  and  final  values  of  the  independent  variable , 
z,  over  which  the  integration  process  is  to  be  carried  out.  The  inte- 
gration step  size  is  also  stipulated.  The  integration  can  be  processed 
in  up  to  3  segments,  each  with  a  different  step  size,  thus: 

ZI     S1  t>  Zl    S2  -Z2     §3  »  zr„ 
or     ZI      S3-  »Z1     S2»-ZF  , 
or     ZI      S1   »ZF  , 

where  ZI  refers  to  the  initial  value  of    z    (in  this  case  equal  to  zero), 
ZF  refers  to  the  final  value  of  z,  Zi  refers  to  intermediate  values  of 
z,  and  Si  are  the  step  sizes.    All  of  these  numbers  are  entered  in 
meters.    The  corresponding  values,  with  decimal  points,  are  placed 
in  the  above  order  in  columns  1-10,   11-20,  21-30,  etc. 
Data  card  #5:    Columns  1-10  contain  the  value  of  the  energy  E  in 
electron-volts,  again,  with  the  decimal  point.    Columns  11-20 
contain  the  quantity  P  which  adjusts  the  magnetic  field  according  to 
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the  base  of  10  amperes  magnetization  current.    For  example,  for  a  field 

of  7  amps,  the  quantity  0.70  would  be  entered  anywhere  in  columns 

11-20. 

Data  card  #6:    Here  are  the  values  of  the  initial  conditions.    The  format 

is  the  same  as  for  cards  4  and  5.    The  first  value  expressed  is  the 

initial  value  of  r  (equal  to  zero),  and  the  next  value  is  the  slope  of 

the  function  at  that  point,  (_^L\    .    This  derivative  may  be  considered 

Vdz/ 

to  be  also  the  tangent  of  the  scattering  angle,  8  . 

Data  card  #7:    This  card  controls  data  print-out.    A  blank  card  will 

suppress  this  output.    Again,  using  the  ten  column  sequence,  this 

format  is  as  follows.    The  first  8  columns  contain  the  title  for  each 

variable  output  desired,  and  the  last  two  columns  contain  the  variable 

identifier.    The  variable  identifiers  are  00  for  z,  01  for  r,  and  02  for 

r'.    An  example:    Z  DIST    00R  DIST    01SLOPE      02 

Data  card  #8:    This  last  data  card  controls  graph  output.    A  blank 

card  will  suppress  this  output.    Remember:    that  only  ONE  graph  for 

each  submitted  program  can  be  drawn  and  it  corresponds  to  the  first 

set  of  data.    The  information  is  contained  in  columns  1-20.    Columns 

1-16  contain  the  graph  title,  and  columns  17-20  indicate  the  variable 

identifiers.    Columns  17-20  must  always  be  0100.    The  origin, 

direction  of  plot,  and  scale  distances  have  already  been  set  within 

the  main  program . 

Running  Time:    The  time  required  for  the  computer  solution  is  dependent 

on  the  step  size  and  the  total  range  of  the  independent  variable.    This 

running  time  has  averaged  about  three  minutes  for  a  single  trajectory. 

7.      Results 

The  outcome  of  this  project  was  the  creation  of  a  single-package 
system  which  generates  hydrogen  ion  trajectories  in  an  axially 
symmetric  field. 

Approximately  100  production  trajectories  were  computed  for 
values  of  energy  E  ranging  from  97  eV  to  450  eV;  magnetizing  current 
from  7  to  15  amperes;  and  scattering  angle  of  35-55  degrees.    The 
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numerical  output  data  and  the  graph  plots  showed  these  orbit  con- 
figurations to  be  quite  smoothly  continuous  and  consistent. 

The  error  in  evaluating  the  intercept  distance  L  was  due  mainly 
in  the  uncertainty  of  reading  the  graph  or  interpolating  the  data  print- 
out.   These  values  of  L  ranged  from  20  centimeters  to  about  80  centi- 
meters, with  an  estimated  error  of  ±1.0%.    The  maximum  radial  distance 
of  the  orbits  was  around  8  cm. 

Comparison  of  these  computer  trajectories  with  those  obtained  by 
the  fine  wire  technique  used  by  Kelly  showed  only  fair  agreement  to 
within  ±10.0%.    The  disparities  in  the  different  sets  of  trajectories 
have  not  as  yet  been  determined. 

This  computer  system  was  checked  with  a  constant  uniform 
(homogeneous)  magnetic  field  against  an  analytic  solution.    The  com- 
parative results  were  excellent  (to  within  ±0.01%). 
Based  on  a  Standard  Trajectory  of     E0  =  300  eV# 

IB  =  Magnet  current  =  12  amperes, 
0   =  45  degrees  : 

-|^-|  T        =    -3.3  ±0.2  cm/deg 

o  8  /  E0 


'  LB 


■^77)  _       =    -9.7    : 0.8  cm/amp 


J  e, 

), 


bE0j  iB,e 


=  0.20  ±0.01  cm/eV 


The  above  derivatives  show  that  as  ©and  Ig  are  increased,  L  is 
decreased.    However,  an  increase  in  E0  produces  an  increase  in  L. 
Each  of  the  partial  derivatives  were  obtained  with  respect  to  a  change 
in  a  single  parameter,  with  the  other  two  parameters  held  constant. 
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APPENDIX  1 
DATA  ON  THE  MAGNETIC  FIELD 

The  magnet  which  provides  this  inhomogeneous  field  was  formerly 
employed  in  the  beta-ray  spectrometer  which  was  reconstructed  here 
at  the  USNPGS  in  the  period  1955-56.    Although  much  was  known  about 
the  magnet  itself,  there  existed  no  data  on  the  magnetic  field  until 
measurements  were  made  by  Kelly  [5  ]  in  1965. 

This  data  has  been  compiled  and  is  herein  available  for  the  field 
strength  both  on-  and  off-axis;  and  also,  for  the  theoretically  computed 
as  well  as  experimentally  measured  values  in  Gauss. 

In  the  diagram  of  the  Magnetic  Field  Strength  in  this  appendix, 
there  are  two  values  listed  for  each  (r,  z) .    The  experimental  values 
in  the  diagram  are  the  upper  ones  listed  in  each  square. 

The  theoretically  computed  values  are  listed  in  squares  and  are 
the  lower  values.  The  method  by  which  these  values  were  computed 
is  as  follows: 

The  field  on-axis  is  given  by: 
H0(z)  =  Bx  +  B2z  +  B3z2  +  B4z3  +  B5z4  +  .  .  .  +  Bl3z12  (3-1) 

(z  expressed  in  meters) 

where: 

B(  1)  =  8.0498907E-02 
B(  2)  =  2.3351831E-02 
B(  3)  =  -4.8664438E+00 
B(  4)  =  4.2887590E+01 
B(  5)  =  -4.3626183E+02 
B(  6)  =  3.9271101E+03 
B(  7)  =  -2.2648265E+04 
B(  8)  =  8.2393030E+04 
B(  9)  =  -1.9359723E+05 
B(10)  =  2.9437706E+05 
B(ll)  =  -2.8025472E+05 
B(12)  =  1.5200646E+05 
B(13)  =  -3.5865075E+04 

Since  A  ^    is  the  only  non-zero  component  of  the  vector  potential,  the 

axial  component  of  the  field  strength,  Hz    ,  can  be  computed  and  compared 

to  values  obtained  by  the  Hall  Probe. 
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v{t(^4 


where    A  ^    is  equation  (5-1). 

And,  after  performing  the  indicated  operations 


Hz    (r,z)  =G1 


ia 


4       -     64  2304 

This  is  demonstrated  in  Appendix  6. 


G4+  . 


(11) 
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APPENDIX  2 

USEFUL  FORMULAE  IN  CYLINDRICAL  COORDINATES 

Conversion  from")  x=r  cos  0 

Rectangular     r  y=r  sin  0 

Coordinates  J  z=z 

The  Gradient 

7 


The  Divergence 

7  ' 


%=1  jT&(rAr)l  +   J  (lA^)+(5Az1 
r    (_    Br      J        r  1 d0    f      \    Sz    J 


The  Curl 


7xA=l[Mi     -   -^rA  0)1 
r  I    d0  dz 


I    Br  BzJ    * 

+     1    [dfrA  0)     _     SAr  ] 

r     [    Sr  B0 


The  Laplacian 

^2 


r\Br    V      dr/J       r2|B02  J      [fcz2    J 
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APPENDIX  3 
POWER  SERIES  EXPANSION  TO  OBTAIN  THE  VECTOR  POTENTIAL 
The  well  known  expansion,  equation  (5-1)  will  now  be  derived. 
This  expression  describes  the  vector  potential  for  an  inhomogeneous 
magnetic  field,  which  acts  only  in  the  azimuthal  direction.     The  field 
on-axis,  H0  (z),  and  its  even  derivatives  can  be  computed. 

Basically,  the  derivation  of  this  expansion  is  as  follows: 
Since  the  field  is  constant  and  involves  no  currents  in  the  region  of 
the  field, 

7  x    B  =  0. 
And, 

"B  =  7*  x    A, 
so  that  the  vector  equation  for  which  a  solution  is  desired  is: 
7   x    (  7   x    A)  =  0 . 
Conditions  for  this  field  are: 

B0    =0, 
and       Ar     =       A  ^    =    0. 


In  cylindrical  coordinates,  7  x    B    becomes 


7  x    B    = 


BB7.       d(rB.) 


a  z 


SB. 


3B 


0> 


OZ 


r     L     ^  r 


and  each  component  must  be  equal  to  zero. 


o  B, 


=     &% 


or 
Br 


oz 


I    [Mz]  _  fa(rA0   )1  =  I    (-  S(rA0  )i 
r     L  B0j       L      Bz        J       r       I         dz         J 


=      dAz 
br 


SAr 

oz 


=    0 


Bz    = 


1     fB(rA^   )1  _    _SAr_    =    1     f    5  (rA  ^      )  1 
r      (_      Br  d0  r      [        ar 


SB, 


(1) 
(2) 


(3) 
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From  equations  1 ,  2  and  3 , 


1-fiA     +    A]         J_[jr     A  _BA_]     =  0 

Br     [  Br  r  J  Bz    j^Bz     r  Bz   J         U' 


and 


and  the  partial  derivative     — —    is  zero  as  the  coordinates  are  indepen- 

o  z 


dent 


Thus 


B2A  I    _BA 

Br* 


A      +     HA 
r     Br  r2  Bz2 


Here,  a  solution  is  assumed  in  the  form  I    an  (z)  rn        where 

n=0     " 

the  coefficients  an  (z)    are  differentiable.    To  determine  the  coefficients, 


BA 


n-1 


f    =     E  nan(z)rn^  , 


and 


and 


and 


-^2-  =    E  n(n-l)an  (z)  r 


n-2 


BA 
Bz 


B2A 


L  an'  (z)rn 


jP"  =    £  V  (z)  rn 

Substituting  in,  and  multiplying  through  by  r2 

£(n)(n-l)an(z)rn    +    In  an(z)rn    -    San(z)rn    +    Ian"(z)rn+2      =0 

E    {  n(n-l)  +  n  -  l}an(z)rn  =    -Ean"(z)rn+2 
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For  n=0#  the  left  hand  side  of  the  last  equation  reduces  to  -a0(z),  and 
there  exists  no  n=l  term.    Therefore, 

-a0(«)  +n|2     (n(n-l)  +n-l}an(z)rn   +nJ^   an"(z)rn+2     =0 

-a0(z)  +    £    {(n+2)(n+l)+n+l}an+2  (z)rn+2     +     £   an"(z)rn+2        =0 
n=0  n=0 

-a0(z)    +    f    {  (n+3)(n+l)}an+2  (z)rn+2    +     L    an"(z)rn+2      =    0 
n=0  n=0 

-a0(z)    +    f    CL(n+3)(n+l)]^+2    (z)    +  an"(z)}rn+2      =    0 
n=0 

The  following  must  then  be  true,  as  for  all  power  series  expansions: 

a0(z)    =  0  ,    hence 

(n+3)(n+l)an+2  (z)    +    an"(z)    =    0 

'    ^+2    (z)    =     -an"(z) 
(n+3)(n+l) 

Whence,  a  recursion  relationship  yields  the  important  needed  para- 
meters.    Obviously,  there  can  be  no  terms  for  n  even. 

For  n  =  1  a3(z)    =    -ai"(z) 

8 

n=  3  a5(z)    =    -a.V'(z)  etc. 

24 

As    A(r,z)    =     £   an(z)  rn 
n=0 

=  a^zjr  +  a3(z)r3  +  a5(z)r5  +  a7(z)r7  +  .  .  . 
Making  the  substitutions, 

A    =  al(z)r-a1"(z)f   +  (-i-J (±)  a, ■■ "  (z)r*  -  . . . 
And,  to  evaluate    ax(z),  solve  for  Bz  using    7  x   A, 
Bz     -^--^    =2al(z)-|V(z)r=   +  ^  a,  ■■  •■  (z)r*  -  . . . 
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Knowing  that  B     is  equal  to  H0(z)  for  r  =  0,  hence, 
H0(z)    =    Bz(r=0)    =    2a1(z)|  .\  a^z)    =    £   H0(z) 

and  the  expression  for  the  Vector  Potential  is  obtained  a; 
A(r,z)    =    -y  H0(z)    - -~-    H0"(z)    + -^    H0""(z)-   ... 
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APPENDIX  4 
THE  DIFFERENTIAL  EQUATION  DEVELOPMENT 

The  essentials  of  this  operation  involve  simply  applying  Newton's 
second  law  in  cylindrical  coordinates  to  the  magnetic  force  acting  on 
the  particle.    Finally,  time  and  the  azimuthal  variables  can  be  elim- 
inated by  means  of  the  two  motion  constants,  energy  and  the  canonically 
conjugate  angular  momentum.     Since  the  particle  will  begin  its  tra- 
jectory on  the  z-axis  ,  the  canonically  conjugate  angular  momentum  is, 
in  fact,  zero. 

Applying  Newton's  2nd  law, 

=    q  (v  x  B)    =    m  (a ) 


mag 


where 


d       r^n 

djT  [v] 


In  cylindrical  coordinates,  the  velocity  vector  is 

_».  .A  .A  .A 

v    =    rr    +    r  0  0  +    zz, 
and,  the  acceleration  vector  is 

a*  =    (r  -  r  02)  r    +    (2r  0  +    r  0)  J  +    (z)  z. 

Knowing  the  vector  potential,  the  magnetic  field    B    and  the 
quantity    v  x  B    can  be  found. 


B    =  7 


«*-7[-w*«)-£tA#3f  +fW 


_d_ 

Bz 


("Is 


But: 


A    =    A0    0 
Ar  =    Az  =    0 


and 

-»- 
B 


-*[-&Hf  +r  -h-(*») 


J3 


Hence: 


* x  f  =  [*  I?  (rA  *)J  *  +  [  7  fr(rA »)  -  7  MrA  *)  * 
+  jj  i_(rA  ^ . 

Equating  each  orthogonal  cylindrical  component,  the  following  is 
obtained. 


m(r  -  r  02)    =    q(0    - — rA  0) 


or 


m(2r  0   +    r  %)    =    q 


[4  U>*)-7  £-M] 


m(z)  =  q(  0   TTrA  0)- 

o  z 


These  three  equations  can  be  reduced  further  to: 


—  (mr)  =  mr  0 2    +    q  0    ^-  (rA  0) 


(4) 

(5) 
(6) 


dt 


(mr2  %  +    qrA  0)    =    0, 


and 


-^r(mz)=qrJ    ^7<A0). 

From  the  middle  equation  above,  it  can  be  shown  that  the  quantity 
conserved  is  the  canonically  conjugate  angular  momentum  about  the 
symmetry  axis  .    Thus, 

Pc       =    mr2  0   +    qrA  0  , 


_       Pc 


qrA 


and, 


F(mr)  " 

* 

Pc 

d    ,    .. 

~t~ 

[  dA 

B2 

*•][>♦  %*] 


[V]  [>-**] 
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Replacing  the  independent  variable  of  time  (  t  )    with  the  distance  on 
the  z-axis  (z)  yields 


and  2       r  t 

(z)2    =     m2      L  qr ^J 

{■  *  m  'J 

This  operation  can  be  performed  as  follows: 

[it  <->]  • 


^mz)    =    z 


v    =    rr    +    r00+zz 


v2    =   t;.  v    =    (f)2    +    (r  0)2    +    (z)2    , 

and 

dr      . 
dz 


h    =    Pr    "     °IrA 


Hence: 

,,2 


Further: 

(z)2    = 


(/) 


(H) 
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and: 


dz2    (z) 


(9) 


'  |_tr<z>J  =    It  117}  |ir  "  AJ  •  (10) 

Combining  the  remaining  equations  8,9,  and  10,  and  then  allowing 
that 


and 


Pc    =    0 
mv 

k  s  T    ■ 

the  following  equation  results: 


L  J   L  J  (4_i: 


where  the  dots  throughout  this  discussion  indicate  derivatives  with 
respect  to  tpne,  and  the  primes  indicate  derivatives  with  respect  to 
the  variable    z.     [4  ] 
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APPENDIX  5 
DEVELOPMENT  OF  THE  VECTOR  POTENTIAL 
To  adequately  describe  the  field  at  every  point,  and  for  use  in 
the  general  trajectory  previously  derived,  equation  (4-1),  the  following 
quantities  are  needed: 

A  =  The  Vector  Potential, 

— —  =    The  first  partial  derivative  of  A  with  respect  to  z, 
o  z 


and 


The  first  partial  derivative  of  A  with  respect  to    r. 


The  vector  potential,  found  by  the  power  series  expansion  in  Appendix 
3,  is 


A  ■  ftGi]-TrtG.]  + aS-cG.]  -7ii^[G*]  +  --- 


r13 

[G7] 


29727129600 

where,  for  the  purposes  of  programming,  the  G  functions  are  defined 
as  follows: 

Gx     =     H0(z) 


G2     = 

d2  tUz) 

Bz2 

G3     = 

S4   H. 

dz4 

G5     ■ 

beHn 

az6 

G        =     *      H' 
G?      ~      S  z12 


V7 


The  first  partial  of   A    with  respect  to    r    becomes, 


br     ~    *[Gi]-  16    [G*]  +   384[G3]"    18432  [G4]+  1474560 


UH rG  -,  .  Ui* 

176047200   L  6J    29727129600 


LGSJ    +     oQvo^nOQcnn      LG7J 


The  first  partial  of   A    with  respect  to    z    is  a  little  more  tedious,  but 
can  be  shown  to  be, 


dA 
bz 


=    T[F>]-    if^^^T^" tF6] 


where  the    F    functions  are  defined  as  follows  (again  for  programming 
convenience) , 


*l  = 

bz 

F2    = 

SG2 

bz 

F3    s 

bG3 

bz 

F4    s' 

bG4 

bz 

F,-     =- 

bGB 

rS 

bz 

F«    =. 

BG6 

bz 

F7    = 

0. 
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APPENDIX  6 
PERTINENT  MACHINE  PROGRAMS 
The  pertinent  computer  programs  used  in  the  solution  of  this 
problem  are  listed  in  this  appendix  in  the  following  manner: 

Program  A    -    Least-squares  polynomial  generation  by  the  PROGRAM 

CURVE,  and  contains  7  pages. 
Program  B    -    Theoretical  off-axis  field  check  labeled  PROGRAM 

WATCHOUT,  and  contains  2  pages. 
Program  C   -    The  main  trajectory  program  labeled  PROGRAM  ORBITS, 

and  contains  15  pages. 
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APPENDIX  7 
GRAPHS  OF  VARIOUS  TRAJECTORIES 
This  appendix  contains  these  four  figures  which  show  trajectories 
as  produced  by  PROGRAM  ORBITS. 


Figure  3.  The  Standard  Trajectory- 
Figure  4  .  Variation  in  Energy 
Figure  5 .  Variation  in  Magnetic  Field 
Figure  6.  Variation  in  Initial  (scattering)  Angle 
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